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Abstract 

We reanalyze the two-loop electroweak hadronic contributions to the muon g — 2 that may be 
enhanced by large logarithms. The present evaluation is improved over those already existing in 
the literature by the implementation of the current algebra Ward identities and the inclusion of 
the correct short-distance QCD behaviour of the relevant hadronic Green's function. 



1 Introduction 



The latest result from the g M — 2 experiment at BNL [Q reported by the E821 collaboration has 
triggered a renewal of interest on the theoretical prediction of the anomalous magnetic moment of the 
muon = 5(5^ — 2) in the Standard Theory Q The attention has been focused mostly on the hadronic 
vacuum polarization contribution, particularly on the accuracy of its determination || [|, ||, g]. 
The major theoretical change, however, comes from a new evaluation of the dominant pion-pole 
contribution to the hadronic light-by-light scattering which has unravelled a sign mistake in 
earlier calculations ^ [0], |l2], [L3| . The correctness of the result in ref. Q was corroborated by 
a renormalization group argument within the low-energy effective theory of the Standard Theory in 
ref. jL4j, with subsequent developments in refs. [jl5|, [l6| and |l7| . The previous calculations have now 
been amended correspondingly |l8|, [ll| ||(| . 

Here we wish to report on a new calculation concerning an interesting class of hadronic contribu- 
tions which appear at the two-loop level in the electroweak sector. They are the ones generated by 
the hadronic 7 — 7 — Z vertex, with one 7 and the Z-boson attached to a muon line, as illustrated by 
the Fcynman diagrams in Fig. 1. 

As first noticed in pl| these contributions are particularly interesting since, in principle, they can 
be enhanced by a large log(M§/mj 2 oop ) factor, where mi oop is the relevant scale in the shaded loop in 
Fig. 1. However, in the absence of the strong interactions, there is an important cancellation between 
lcptons and quarks within a given family, as a consequence of the anomaly-free charge assignment in 
the Standard Theory |^2|, . The purpose of the present work is to analyze what is left out of this 
cancellation in the sector of the u, d and s quarks, where the strong interactions turn out to play a 
subtle role. 

When dealing with the strong interaction effects, it has become a common simplifying practice 
among some theorists, to assume that the main effect of the strong interactions is dual to a modification 
of the u, d, s quark masses in the QCD Lagrangian, to "constituent-like" quark masses of the order of 
300 MeV for the u, d quarks and 500 MeV for the s quark. In this way the contributions from Fig. 1 
were found to be rather small j22|, |23| 0. However, in view of the expected accuracy which the new 
BNL experiment will eventually reach, it becomes necessary to have a more reliable determination 
of the size of these contributions. This has now become possible within the framework of the 1/N C 
expansion in QCD. There are indeed recent theoretical developments of this non-perturbative analytic 
approach, which have been applied mostly to the calculation of hadronic weak processes 0, but which 
turn out to be useful as well for the evaluation of this particular class of hadronic contributions to a M . 
This paper is dedicated to the evaluation of these contributions. 




Fig. 1 Fcynman diagrams with the hadronic 7 — 7 — Z vertex which contributes to the muon anomaly. 
1 For a recent review see e.g. ref. |s| where some of the earlier references can also be found. 

2 Notice that ref. |22| us es a more sophisticated version of this approach where at least chiral symmetry breaking is 
correctly implemented fajj. 

3 See ref. 1251 for a recent review, and refs. Pq| to |30| for details. 
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The relevant terms in the Lagrangian of the Standard Theory which we shall be needing are the 
following: 

la. The Leptonic Neutral Current Lagrangian 

/Heptojis _ Tom Af* £ t(0) 7V 

^no 2 cos 9\\? ' 

where 

/ / 

with 

v/ = Tsf — 2Qf sin 2 9w , and a/ = T3/ . 
In our case, we shall only need the muon components 

Qn = -1, and a /i = -l/2; 

i.e., the muon electromagnetic coupling and the axial coupling of the Z to the muon line: 

-ep,(x)j a fj,(x)A a (x) and - — -fi(x)j a j 5 n(x)Z a (x) . 

Z COS Uw 2 

lb. The Hadronic Neutral Current Lagrangian in the sector of light quarks 

1 - 75 , (0) 1 + 75 



where q = (u, d, s) and 



,(0)i^_/5, (0). 

" 2 " 2 



q(x) . 



i<°> = eQ L [Aft - tan 8 W Z^ + \ Qfz» and rjW = eQ fl L4 M - tan 9 W Z^] . 

Z COS u w 

In our calculation, only the hadronic electromagnetic coupling: 

e g(x) 7 ^gg(x)A A1 , Q = Ql=Qr = diag (2/3,-1/3,-1/3) ; 
and the hadronic axial coupling: 

~\ o \ q{x) 1 ^ 5 Qfq{x)Z ll {x) , Qf = diag (1, -1, -1) , 

Z Z COS uy/ 

will be needed. 

1c. The coupling of the unphysical neutral Higgs field to the light quarks and to the muon 

Uo = T -? ir <S> (x)J* (x), (1.1) 
2 cos uw 

with (a„ = 1/2 , a d = a s = -1/2) 

J®{x) = Z 22 a qj^Q( x )h5q(x) + 2a, fi j^fl(x)ij 5 n(x) . (1.2) 

g— u,d,s 

The hadronic electroweak vertex which we have to compute, in the notation corresponding to Fig. 1, 
is the following Q 

W)\VrW\l*(P)) =u(p>)r p (p>,p)u(p) = He)(-ie) fj-^z-) fj-^-) X 
p \4cosfljy/ \4cos6'vk/ 



4 We use the following conventions for Dirac's 7— matrices: {7^,71/} = 2g M „, with g M „ the Minkowski metric tensor 
of signature {+, — , — , — }, 75 = «7 7 1 7 2 7 3 , and the totally antisymmetric tensor e M i, pCT is such that eoi23 = +1 ■ 



(2tt) 4 g 2 (p'-p-g) 2 -Ml 



4-m, 



-7 75 +7 75- 



u(p)x 



dfixe*-* / d?yjV^^n\T{Vr(x)A?(y)Vr(0)}\n)> 



(1.3) 



where 



^ / ! m (a ; ) = <?W7f* Q 9(2;) > = q(y)iui5 Q L i(y) 



and |fl) denotes the full QCD physical vacuum. Here, for the purpose of simplicity, both the photon 
and Z propagators in the second line of Eq. (1.3) are evaluated in the Feynman gauge. We dedicate, 
however, Appendix A to a discussion of questions related to gauge dependence. The anomalous 
magnetic moment contribution from this vertex is then defined by the corresponding Pauli form 
factor at zero momentum transfer, i.e., 



F 2 (0) = lim tr{{f + m ^)k p 2 {p',p){f + m^)T p {p' 1 p)} 

k 2 —>0 



(1.4) 



where (p' — p + k) 



k 2 Am 2 — k 2 



m„ 1m„ + k 
vP _rrv j, ; (p +p 'y 



k 2 (4m 2 - fc 2 ) 2 



(1.5) 



is the projector on the Pauli form factor, and m p denotes the muon mass. 



2 The Master Green's Function 

The contribution we want to compute is governed by the hadronic Green's function 
W^ p (q,k) = [ #xe** I d i ye^^y(n\T{V-(x)A-(y)V-(0)m, 



(2.1) 



with k the incoming photon four-momentum associated with the classical external magnetic field, as 
illustrated in Fig. 2. 




Fig. 2 Symbolic representation of the (VAV) three-point function in Eq. (2.1) 



The general Ward identities which constrain this type of three-point functions in QCD are discussed 
in Appendix A. In particular, W pvp {q,k) has an electroweak U(l) anomalous term. It is this term 
which, when lumped together with the corresponding contributions from the c quark and the e and fi 
leptons to the 7 — 7 — Z loop, gives a finite gauge invariant contribution to a M which was calculated, 
independently, in refs. |2^, |2^|. Here we are particularly interested in the contribution to a p from the 
non-anomalous part of W pvp {q, k), denoted by W liVp {q, k), i.e. 



W^ vp (q, k) = W livp {q, fc) anomaly + W liup {q, k) . 



(2.2) 
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where 



ah 



12^ 3 (g - A;) 2 ^ ' 



(2.3) 



The function W p „ p (q, k) obeys the trivial Ward identities 



q fi W^ p (q,k) = ) and k" W^p (q, k) = , 



(2.4) 



but, as discussed in Appendix A, it has both a longitudinal and transverse component in the axial 
neutral-current i/-index: 

W^ p (q, k) = W^(q, k) + W^(q, k) . 



Differentiating the second identity in Eq. ( 2A ) with respect to k we get 



W^ p (q,k) = -k^W^(q,k). 



(2.5) 



As we shall see in the next section, extracting one power of k from W^ p s (q, k) is all that is needed in 
order to obtain the corresponding contribution to the Pauli form factor at zero momentum transfer, 
i.e. the anomalous magnetic moment a p . Therefore, without loss of generality, we can set 



WX°(q, k) = k°W pvpa {q) + 0{k 2 ) 



where 



W pupa {q) = 



dkP 



W pva {q, k) - /^_Sa (q - k) v 'W^ a (q, k) 



(2.6) 



The most general pseudo-tensor W pupa (q) which satisfies the Ward identities is then constrained to 
have the form (Q 2 = —q 2 ) 



( q ) = i W(Q 2 )[q 



Qa^piyapQ 



(2.7) 



As we shall see in the next section, the contribution we are looking for can be expressed as a simple 
integral of the function W{Q 2 ) over the range < Q 2 < oo in the euclidean Q 2 variable. We shall 
show in section 4 that, in the chiral limit, this function falls as 1/Q e at large Q 2 and, in the presence 
of massive light quarks, the large-Q 2 fall-off of W(Q 2 ) goes as l/Q 4 only. 



3 The Contribution to the Muon Anomalous Magnetic Moment 



Since what we want is ^(O) we need only to keep the projector in Eq. ( jl.5| ) to its simplest form in 
an expansion in powers of momentum fc, 

A 2(P'P) = 4^2 h 



2m,, 



This is because the product 

(p + m p )A p 2 {p',p)(p" + m p ) 



1 

4fc2 ( 



: + m ll ) 



-ft- 



k" 



p' 



O(k 



5 There is a priori another possible tensor structure which also fulfills these requirements: 
However, using the identity 



£ p,vpa \ _|_ £ vptr\ n _|_ e pcr\p, u _|_ e <y\^v ^p £ \pv p _ q 



it can easily be shown (contracting the previous identity with qx) that it is identical to the one in Eq. (2.7). 
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which then appears in Eq. (1.4), is proportional to one inverse power of k . On the other hand, 
gauge invariance forces the Green's function W pl/p (q, k) to be proportional to one power of k at least. 
Altogether, this gives the required powers of k which are needed to define ^(0). We can, therefore, 
simplify significantly the calculation of Eq. (1.4) to that of the simpler expression: 



F 2 (0) = (-e 2 ) 



lim 



d 4 q 1 



16 cos 2 9 W Mf k^a J (2tt) 4 q 2 \q 2 - M. 



M 2 



4fc 2 



— 5 tr Ut + m,>) 



Y ft - [ k» + ft 

m,. 



{i> - 4+ m M ) v v ( 
7 o ^ —1 75 + 7 75- 



+ 



N c 4q u 



12^ 2 3 g 2 WQ 



q 2 — 2q-p q 2 + 2q-p 

a q a k a - iW x ;«l (q, k) + k a [q p e pvaa q a - q a t pvap q a \ W(Q 2 ) 



(3.1) 



where the first term in the last line of Eq. (3.1) is the contribution from the anomaly to the (VAV) 
Green's function, already discussed in ref. |22| . The second and third terms are the non-anomalous 
contributions from the longitudinal and transverse components of the same Green's function. Let us 
first simplify, still further, the last contribution. 

The integral over the four-momentum q from the last term (the transverse component), is conver- 
gent in the infrared because W(Q 2 ) has no Goldstone pole at Q 2 — > and in the ultraviolet because 
of the QCD short-distance behaviour of W(Q 2 ), which we shall later discuss. Therefore, to leading 
order in powers of the lepton mass, we can perform the trace and then the integration over the angles, 
setting 

1 1 

?> (3.2) 



q 2 ± 2q-p 

in Eq. (3T). The integral over the solid angle dQ q becomes trivial, and the transverse contribution is 
then given by the expression 



F 2 (0)1 
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lim 



16 cos 2 W M 2 fc2->o 16?r 2 



dQ z Q 



Q 4 \Q 2 



Ml 



Q 2 W(Q 2 )x 



-1 
8k 2 ' 



tr < (jS + m p ) 



7 P # 




[7^ + m p )j"j 5 + 7 I/ 75(^ + "fy)7 M 



M 2 Z 10g m 2 



The remaining algebra leads to a remarkably simple expression: 

,2 



*a(°)L 




dQ 2 



M 2 Z 



Q 2 + Mj 



W(Q 2 ). 



(3.3) 



(3.4) 



The integration of the anomalous term in Eq. (3.1) is in fact slightly more complicated, because 
there we cannot simplify the muon propagators as in Eq. (3.2). It is then useful to combine denomi- 
nators using two Feynman parameters: 



{q 2 ) 2 (q 2 -M 2 z ){q 2 ±2q-p) 



= 6 



dxx 



dyy- 



1 



where 



[{q-lf-B?T 
- Mf x{l-y)+m 2 (l-xf . 



Then, after the shift q 

m ( m u i Af 



I = ±p(l - x) and R z 

q + I of the integration four-momentum, and upon neglecting terms of 



one obtains the result (in the Feynman gauge) 



■5 



dx x 2 



dyy 



1 + 3(1 



(1 



x(l-y) + ^-(l-z) 



M 



2 r 



\/2 87r 2 7r c 



4 , M i 



G F a 

y/2 8TT 2 7T 



x 18.69. 



^l-^ + ^l-x) 







A]| l0g 



Ml 



(3.5) 



in agreement with the result |^| in ref. ^2|. The evaluation of F^O) in an arbitrary gauge, with a 
discussion of the physics behind it, can be found in Appendix A. 



4 Constraints on the Function W(Q 2 ) from the Operator Product Expansion (OPE) 



In order to perform the integral over Q in Eq. (3.4) we require information on the function W(Q ). 



We shall now discuss the constraints on this function, at large Q 2 values, which follow from QCD 
at short-distances. For pedagogical reasons, we shall first discuss the hadronic contributions in the 
chiral limit, where the light quark masses are neglected. However, the leading effect due to the explicit 
chiral symmetry breaking induced by the finite light quark masses, turns out to be important and 
rather subtle. It is at the origin of a conceptual discrepancy with the treatment reported in ref. ]23[ ]. 
We shall, therefore, discuss as well the contributions from explicit breaking, later, in section 6. 

Since one of the photons is on shell, the relevant short-distance behaviour we are looking for is the 
one giv en b y the OPE of the two currents V° la (x) and A'^(y) in the three-point function W tlup (q 1 k) 
in Eq. (|2.l|) when x — > y, because these are the currents where the hard virtual momentum q can get 
through |3lf| ( from here onwards, the subscript non-anomalous in T-products will be understood). 
This three-point function is represented symbolically in Fig. 2. The quantity we are interested in is 
defined by the limit 

ZV(<?)= lim / d 4 ze^ 2 T{l/;'"(z)^(0)} , (4.1) 
and the relevant term in the OPE of these two currents is the one with a tensor structure 

, 0^(0 ) 
(Q 2 ) 



Uy.v{q) = i [qst^ a pq a - q^^ va &q a \ C /nt2 . p H — , (4.2) 



where C is a dimensionless constant to be calculated, p denotes the lowest possible power and the 
ellipsis stands for other contributions proportional to higher powers of 1/Q 2 and/or to other tensor 
structures which cannot contribute to the muon anomaly. The dimension of the 0@ s operator is such 
that: 

d [O 136 ] -2p = 0; 

i.e., 0@ s must be an operator of even dimension, and antisymmetric in (38. 



There are in fact two types of potential contributions to the r.h.s. of Eq. (4.2) which, physics- 
wise, are generated by the emission of soft virtual gluons and hard virtual gluons in the T-product in 
Eq. (f4.ip. As discussed in Appendix B, there is no possible contribution from soft gluons, on symmetry 



grounds, with the tensor structure of Eq. (4.2) and with a power p < 3. The leading contribution, 
which appears with a power p — 3, comes from the perturbation theory expansion in hard virtual 
gluons and it is the one shown in the diagrams of Fig. 3 below. 

6 The constant term inside the bracket in Eq. (20) of ref. |22| should be 2/3 instead of 4/9. 
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1 Mf" 

(2) 



V 



Y \\f' 




Fig. 3 Feynman diagrams showing the hard gluon exchanges which contribute to the leading 1/Q e 
behaviour of the operator (q) in Eg. . 



The details of the calculation of the contribution from hard virtual gluons to U lll/ (q) in Eq. (4.1) are 
explained in Appendix B where it is shown that, 



iv {q) = i [gse^ a pq a - qpt^asq*] (-27r 2 — ) 



0^(0) 
(Q 2 ) 3 



(4.3) 



and where the operator O^ 5 (0) is found to be the four-quark operator 

~2 



0^(0) 



o x u(T^u) (uu) + \ (da 0S d) (dd) + \ (sa^s) {ss) 
o o o 



(0). 



We find, therefore, that the constant C in Eq (4J2), to lowest order in perturbative QCD is 



C 



7T 



We are now in the position to evaluate the high-Q 2 behaviour of the invariant function W(Q 2 ) in 
Eq. (2.7). For that we first insert the result obtained in Eq. (4.3) into the r.h.s. of Eq. (pi|). What 
appears then is the two-point function 

¥ 3S p (k) = f d^ye- lk -y{n\T{0^{Q)V-{y)} |fi) . 



In our case, the four-momentum k in this two-point function is a soft momentum] and in fact, what 
we need for our purposes is only the term linear in this soft momentum, i.e. 



d 
dk~P 



fe^O 



i / dSv P (n\T{o^(o)vr(y)}m 



(4.4) 



*$V)„ f d 4 y y p (fi \T [Vr(y) (q(Q)o- f3S QQfq{0)) } |fi) , (4.5) 



where in going from Eq. (4.4) to Eq. (4.5) we have used the l/iV c -expansion, keeping only the leading 
contribution, and SU(3) symmetry for the single flavour vacuum condensate (tptp)o- Our problem is 
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then reduced to the one of evaluating the behaviour of the (VT) correlation function which appears 
in the r.h.s. of Eq. ( ff.5| ) at a small momentum transfer. 

On general grounds (conservation of the vector current, parity invariance and SU(3) symmetry) 
and in the chiral limit, the generic (VT) correlation function 

d 4 y e* k -y(n |T | UlJ^A (V) (V^y^) (0)} = [hPsi - k 5 5?) S ab U VT (k 2 ) , (4.6) 

depends on one invariant function HvT(k 2 ) , where A" denotes the Gell-Mann matrices, normalized 
as tr(A a A b ) = 2S ab ', with A = (^2/3) 1 , and a = 0, 1, ... 8. Using the relations 



A 3 1 A 8 
Q = — + ^—, and QQ^> = 



(3> = I(*! + J_*! 
2 ' 2 ' ^ L 3V2 V32 ' "V 3 2 



we can write our result in Eq. fl4.5| ) in terms of Hvt(0) only, as follows: 

"^^Lo = l{ 1 + l) (WhW - ^^)nvr(o) . 

The function UyT(k ) , in the minimal hadronic approximation (MHA) to large- N c QCD, has 
been evaluated in ref. [p2|. In this approximation, it has the simple pole form 



n ^(fc 2 )| MHA =-W) 



1 



fc 2 - M 2 



with M p the mass of the lowest vector state in the large-7V c QCD spectrum. The asymptotic behaviour 
of the function W(Q 2 ) at large Q 2 values is then fixed, in the MHA, with the result 



lim W(Q 2 )\ 



16 2 a s 1 (ip\j})l 
9 ^ ir Ml Q 6 



(4.7) 



5 The Function W(Q 2 ) in the MHA to Large N c QCD and in the Chiral Limit 

Since W(Q 2 ) has dimensions of an inverse squared mass, it is convenient to define 

1 



W(Q 2 ) = ^w(s). with : 



Ml 



(5.1) 



The function w(z) in large— iV c QCD is a meromorphic function in the complex z-variable. It has the 
pole structure shown by the Feynman diagrams in Fig. 4, where the denote vector narrow states 
and the axial-vector states. 



V 



V 



R 



v 

x A V 

— X + X" 



V 

X 



R 



V 

A v x 
■x -|- x — 



'J 



J 



Fig. 4 Feynman diagrams which contribute to the function w(z) in Eq. (5.1) in large-N c QCD 
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Using partial fractions, we can write in full generality 

. . \ -> ctR . Mn 

and Mr the mass of the R narrow state. The short-distance constraints on the unknown residues an 
are 

ctR . — , a/?ZR = , and ^ a_R2^ = TZ , 

R R R 

with TZ fixed by the leading term in the OPE expansion 

lim w(z) = Tz\ +0(4r)- 

2— >oo 2* 

Therefore, in large-iV c QCD 

The minimum number p of poles required to have the same asymptotic behaviour as the OPE 
result in Eq. (4.7) is p = 3. This minimum number of poles corresponds in this case to the minimal 
hadronic approximation to large- N c QCD, which we have used already in several other calculations Q. 
In this approximation, the expansion of the meromorphic function w(z) in rational functions ^] is 
limited to the first three poles, and can be written as follows: 

, , ai Pi a-2 „s 

w( z ) = — T + i + : 7 » (5-3) 
z + l z + — z + r 

9A 

with ai, Pi and a 2 unknown residues, which we shall determine readily, and 

g A = M 2 p /M 2 Al and r' = M 2 p ,/M 2 p . 
As we have seen, the OPE evaluated in the MHA, fixes the value of TZ to 

W= "^2t<M. (5.4) 

There follow then three short-distance constraints on the three unknown residues cti, a 2 and p\: 

ot\ + Pi + a 2 = , 

oti H Pi + r'a 2 = , 

9A 

1 , 2 

Q'i + -5- pi + r a 2 = TZ , 



with the solution: 



a 2 = -72. 



(l-. 9A )(r'-l) ' 
9A 



(l-g A r')(r'-l) 



Pi = TZ- ^ A 



(1 - 9a)(1 - g A r') 

7 See ref. |25| for a recent review, and refs. |26| to po| for details. 

8 Recall theMittag-Leffler theorem for meromorphic functions, see e.g. ref. section 7.4 
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In particular 

w(0) = a x + gA0i + —^i = Tl-^ 



The function w(z) is then fully determined and the integral which defines ^2(0) in Eq. ( p.3| ) can be 
done. The shape of the function w(z) normalized to its value at the origin is the continuous red curve 
shown in Fig. 5 below. As expected, it is a sharply decreasing function in z. 

At this stage, it is perhaps useful to compare the function W(Q 2 ) we find in the MHA to large--/V c 
QCD with the corresponding function, evaluated in the constituent chiral quark model (%QM), which 
was used in ref. |22| to obtain ^(0)^ as an estimate of the error in ^(O). The function W(Q 2 ) in 
the xQM has a simple parametric form: 

W X QM(Q ) = ^9A J q d X J Q dy M , +Q2y{i _ y) , (5.5) 

where gA — 1/2 is the axial-coupling of the constituent quark^. With the same normalization as in 



Eq. (y), i.e., 

1 Q 2 
W x q M (Q 2 ) = j^ w xQm( z ) > and z = > ( 5 - 6 ) 



we have 



Nr. 



I Ml 

12tt^3 9A M 2 X 



*8 



M 2 



tlx 



1 x , xy-y(l- y) 
dy- 



Q 2 



121 



?y(i - v) 



Tog 



VMf. 



iM- 



(5.7) 



The asymptotic behaviours of this function are: 



— N M 2 




and 



lim w x q M (z) 



12tt 2 3 j Q 2 



. Q 2 I 8 Ml 



The xQM not only fails to reproduce the large Q 2 behaviour of the OPE in QCD, but it has the 
opposite sign to the MHA to large-A^ c QCD in the chiral limit. Notice that, even though one is 
interested in the large-Q 2 limit, there is one photon whose momentum k 2 vanishes (see Fig. 2). This 
makes the Green's function W(Q 2 ) in Eq. ( |2.7| ) a nonperturbative object which cannot be calculated 
in terms of free quarks, such as those of the xQM. 

The shape of the function w^Qy(z) normalized to its value at the origin is the dashed curve plotted 
in Fig. 5. It is yet another example of a constituent chiral quark model prediction which deviates 
substantially from the short-distance behaviour of QCD. The resulting contribution to i<2(0)| in 
the xQM is as follows: 



9 See the second reference in 



G 

V2 Stt 

G F m 



ml a 



Nr. 



7T 3 



2 M§ 4 

3 S MK ~ 3 



O 





, M 2 z\ 








log 4J. 



x 3.08 , for M Q = 330 MeV and g A = l/2. 



(5.8) 
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Amusingly, the integral in Eq. (5.2), in the MHA to large- N c QCD can be easily done analytically. 

, ' M 2 

Upon neglecting terms of O 



Mi 



the result reads 



y/2 2^ 



K 



9a 



{l-g A r')(r' -1) 



logr' 



9\ 



(1 - g A )(l - g A r') 



Numerically, for g A = 1/2 and r' = 3.5, we find 

^(°)iir=- 



V2 2^ 



K x 0.128. 



(5.9) 



In order to evaluate, in terms of hadronic parameters, the short-distance residue 1Z in Eq. (5.4), we 
observe that the same factor ~ 2a (^V)o appears in the OPE of the (LR) correlation function where 
it has been shown |M that in large-A^ c QCD 



4tt 



ml 



ii 



In this case, the MHA corresponds to a saturation with a vector state and an axial-vector state only, 
with the result 



-in 



1 



167T 



(5.10) 



Using this relation, as well as the familiar vector dominance expressions |3E 



/* = 2 ^, & = and M$ = 9aM\ 1 , 

1V1 A 1 



Mj 



which as explained in ref. |2j| can also be viewed as predictions of the MHA to large-iV c QCD, our 
hadronic estimate for the residue TZ is then 

« 8 *o 
7£ = — 

9 M 2 p ' 

with Fo the pion coupling constant in the chiral limit (Fq — 87 MeV). This estimate allows for a neat 
comparison with the contributions calculated in refs. j22l E3[, since we now have 



F 2 (0)| 



G F m 2 . a 
V2 8tt 2 tt 



(5.11) 



Even allowing for a substantial error in our determination of (ipip}o hi Eq. ( 5.10| ), w e ob tain, in the 
chiral limit, a much smaller contribution than the one predicted by the xQM m Eq. (5.1 



6 Beyond the Chiral Limit 



As explained in Appendix A, the axial Ward identity in Eq. (A.3), in the presence of light quark 
mass terms in the QCD Lagrangian, is no longer given by the anomaly alone. As a result, the non- 
anomalous part of the Green's function W llup {q,k) in Eq. (2.1), which we have been denoting by 
W tivp {q,k) (see Eq. (|2.2[)) is no longer fully transverse. It acquires a longitudinal component 



where 



(6.1) 



d 4 xe^ x J d 4 ye- ik - v (n\T{V™(x)V™(y)J* (0)} | Q) , 



(6.2) 



and J* is the current defi ned i n Eq. (1.2), where the quark masses enter, linearly, as explicit couplings. 
It is the third term in Eq. (A.2) that we shall now be concerned with. This term, when added together 
with the contribution from the unphysical Higgs coupling in Eqs. (1.1) and (1.2) produces an extra 
contribution to the muon anomaly 



F 2 (0)1 



9 



lim 



16 cos 2 W M 2 Z ki-,0 J (2tt) 4 q 



d 4 q 1 



2fc 2 



tr 



P 
m 



q 2 - 2q-p 



-75 +75" 



q 2 + 2q-p 
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2M 



< 

q— u,d,s 



d 4 ye- lk -y(n \T{VT(x)V?»(y) g(0)z 75 g(0)} | 0} 



(6.3) 



which we are next going to evaluate. 

At the level of accuracy that we are interested in, it is su fficie nt to consider the contribution from 
the strange quark in the sum of light quark flavours in Eq. (6.3), and let m u = vtid — > . Then, the 



relevant hadronic three-point function, in the notation of ref . 1 32 



is 



(IW^fe-fe) = J d*xe^J d 4 ye-' lk -y(n\T{V-(x)V-(y) s(0)z 7sS (0)} | f7> 

= ^pafUfi-tf) (fj Uv (q 2 , k\ (q - k) 2 ) . (6.4) 



In the MHA to large-A^ c QCD, it is found [EM M that: 



HT A (q 2 ,k 2 ,(q-k) 2 ) 



q 2 + k 2 + (q-k) 2 -^ 



(q 2 - M 2 )(k 2 - M 2 )[(q - k) 2 - M 2 } 



(6.5) 



where, in the denominator, we have kept a pseudoscalar Goldstone mass M 2 in order to handle the 
infrared logarithmic divergence which, otherwise, would appear after integration over the ^-momentum 



in (6.3). Here, we are concerned with the limit of Eq. (6.5) when k is soft. Then 

(#)o 



lim^r (q 2 ,k 2 ,(q-k) 2 ) = 



N c Ml 



N c My 



M 
Mr. 



fe^O 



M 2 



M 2 



M 2 



M 2 



(6.6) 



The integral over the momentum q can now be made, very much the same way as the calculation of 
the anomalous contribution in Eq. (3.5), with the result 



G F mla 4 167r 2 m s (<^)o 
F 2 (0)|, ong = - — ^- -- 



1 



M 2 
Ml 



1 r l 

2 



dxx 



dyyx 



Nc M 2 

87r2 F i J x{l - y) + ^-(1 - x) 2 



(1-x) 2 



My 



N c M 2 



M 



f2 



M 2 



, 8^ F 2 M 2 J M 2 x (i _ y) + g(i _ x ) 



1 - -r 



vx(l-y) + ^(l-xf 
(1-x) 2 



M2 x(l-y) + ^(l-x) 2 



(6.7) 



In the approximation where m 2 -C M 2 , to 2 <C M 2 , and keeping only logarithmic terms and 
constant terms, we get 

G F m 2 a 4 8vr 2 F 2 



-^2 (0)| long - — 7=-^~2~ n 



y/2 8tt 2 tt 9 M 2 



v 



N c My M 2 K ( M 2 1 
8^ ~Ff m 2 " I ° g ^| + 2 / 



iV c My 



8tt 2 i^ 2 



M 2 

K 



M>M*-M* 1 ° S W\> (6 - 8) 



M 2 



where we have used the current algebra relation m s (ipip)o — ~FqM^ . For the numerical evaluation, 
we take F = 0.087 GeV, M K = 0.498 GeV, M v = 0.770 GeV and let M vary in the range \M 2 K < 
M 2 < wMft, which includes the values of the 77 and rf masses induced by m s ^ 0. This results in the 
value 



(6.9) 
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where the first error is the one from the variation of M and the second one is a 30% systematic error 
from the MHA approach. This result represents a positive contribution, smaller than the contribution 
from the anomaly, but much larger than the one from the transverse part of the hadronic (VAV) 
Green's function in the chiral limit. 

Finally, we still have to consider possible chiral corrections in the evaluation of the transverse 
term W^p 3 (q, fc); i.e., terms linear in quark masses which may contribute to the invariant function 
W(Q 2 ) in Eq. ( |3.4[ ). The chiral corrections which are potentially important are those which modify the 



asymptotic behaviour of the function W{Q 2 ) in Eq. (4.7) to a less sharp behaviour. The calculation 
of the leading behaviour, 1/Q 4 in this case, is technically rather involved and we have relegated it to 
the Appendix C, where it is shown that 

Hm m o 2 ) = --— (4w " = md - m ^)o 16 2 as j_ mi ( 

q™=o WJ 9 A/ 2 Q 4 + 9 ttA/2 Q6 > ^- w ) 

where chiral corrections to the 0(1/Q 6 ) term and higher order terms have been neglected. 

With W(Q 2 ) = -p^-w(z) and z = , we now have to construct the corresponding MHA to 

large-7V c QCD for the meromorphic function w{z). From the leading OPE behaviour above, it follows 
that the minimum number of poles required is now p = 2. The resulting w(z) function is then 

/ \i Km .., m 4 (4m u - m d - m s )($ip) 
w(z)\ = - — r— , with 7Z m = — — ——. , (6.11) 

which leads to an anomalous magnetic moment contribution 

„, n ,,MHA G F m 2 a g A . ( 1 \ ^ 
itra M ^ 2nl-g A \g A J 

Neglecting m u and with respect to m s , and using again the current algebra relation — m s (i/''0)o — 
F%M\ , we obtain (with g A = l/2, F = 0.087 GeV, M K = 0.498 GeV) 

F 2 (0)r A c%4- log2^#^| = %4-x 0.065. (6.12) 



In this case, however, we know that the MHA in Eq. ( 6.11 ) fails to reproduce the 0(1/Q ) term 
in the OPE, which we know is there, even in the chiral limit. In order to take into account this 
information, we can go now one step further beyond the MHA and incorporate three poles, like in 
Eq. (^^) with the residues constrained now by short-distance behaviour as follows: 

+ Pi + 0-2 = , 

ati + — fix + r'a 2 = -TZ m , 

gA 

l 2 

ai H — 5-P1 + r a 2 = TZ . 
9 a 

The resulting w(z) function in this improved case is shown in Fig. 6, in 10 -3 units (the continuous 
curve). For the sake of compariso n, we also show in the same figure and in the same units the curve 



corresponding to the MHA in Eq. ( 6.11 ) (the dashed curve). The integral over z in the improved case 



can still be done analytically and, neglecting terms of O {^pr log |^fj , we get 

, ,|HA3- Pn ies = G F ml g f [n m + (n + K m )g A ] log/ [TZ+(l + r')1l m }g 2 A logg A 



V2 2tt t (l-g A r'){r'-l) (1 - g A ){l - g A r') 

Ge- rn 2 a 

-4 - x (0-04 ±0.02). (6.13) 
72 8?r 2 tt 
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As expected, the effect of including the 1/Q 6 terms in the asymptotic behaviour of the W(Q 2 ) function 
reduces the integral, but the sign is still the same as in the MHA approximation result in Eq. (6.12). 
The final error we quote is a very generous estimate of the systematic errors. 




7 Results and Conclusions 

The normalization reference of the electroweak contributions to the muon anomaly is the one-loop 
result M, M, M M, M 



V2 8tt 2 



1 



1 - 4suT 



V) + Olj£ log 



Ml 



O 



Ml g ml 



where the weak mixing angle is defined by sin 2 9w — /Mz ■ Numerically, with G F = 1.16639(1) x 

1CT 5 GeV -2 and sin 2 6 W = 0.224, 



19.48 x 1(T 10 . 
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while we recall that the present world average experimental error in the determination of the muon 
anomaly is 

Aa u \ = ±15.1 x 1CT 10 , 

J 4 I Exp. ' 

and, hoping for a continuation of the BNL experiment, it is expected to be further reduced by still 
a factor of four. A theoretical effort on the evaluation of the two-loop electroweak corrections is 
therefore justified. 

It is convenient to separate the two-loop electroweak contributions into two sets of Feynman 
graphs: those which contain closed fermion loops, which we denote by a/f w ^ (ferm) , and the others 
which we denote by a^ W ^ 2 ' (bos). In this notation, the electroweak contribution to the muon anomalous 
magnetic moment is 

of = + af^bc.) + <^)(fe rm ) . (7.1) 

We shall review the calculation of these two-loop contributions separately. 



7.1 Bosonic Contributions 



The leading logarithmic terms of the two-loop electroweak bosonic corrections have been extracted 
using asymptotic expansion techniques (see e.g. ref.p[.) In fact, these contributions have now been 
evaluated analytically, in a systematic expansion in powers of sin 2 9w, up to 0[(sin 2 #vk) 3 ] , where 

log ^72- terms, log -rjf- terms, -jft log jjt- terms, terms and constant terms are kept p3|. Using 

W H W H 

sin 2 W = 0.224 and M H = 250 GeV , the authors of ref. Gl find 



EW(2) 



(bos) = ^ ^ 2 



-5.96 log 



M 2 
lvl w 



0.19 



a 

72 8tt 



m 2 . a 



= ^-JL-x (-78.9). 



(7.2) 



7.2 Fermionic Contributions 

The discussion of the two-loop electroweak fermionic corrections is more delicate. As already men- 
tioned in the introduction, because of the cancellation between lepton loops and quark loops in the 
electroweak U(l) anomaly, one cannot separate hadronic effects from leptonic effects any longer. In 
fact, as discussed in refs. |^2[ [2^|, it is this cancellation which eliminates some of the large logarithms 
which, were incorrectly kept in ref. It is therefore appropriate to separate the two-loop elec- 

troweak fermionic corrections into two classes: One is the class arising from Feynman diagrams like in 
Fig. 1, with both leptons and quarks in the shaded blob, including the graphs where the Z lines are 
replaced by <i> lines, if the calculation is done in the £z-gauge. We denote this class by a^ W< " 2 \l, q) . 
The other class is defined by the rest of the diagrams, where quark loops and lepton loops can be 
treated separately, which we call it a^ W ^ (ferm-rest) i.e., 

a ^( 2 )(fcr) - a^ 2 \l,q) + a* W ^ (ferm-rest) . 

The contribution from a^ W ^ (ferm-rest) brings in to 2 /M^ factors. It has been estimated, to a very 
good approximation, in ref. p3j with the result 



EW{2) (f ,\ G F "V a v 

a,, I term-rest I = — = — — - — X 

M y/2 8?r 2 tt 



5 m? , to? 7\ 



2sin 2 6W V 8M 2 , ™° 3 



where A H ig gEi denotes the contribution from diagrams with Higgs lines, which the authors of ref. 
estimate to be 

A HlKgs = -5.5 ±3.7, 



and therefore 



af L W ^ (ferm-rest) = ^ ^ - X (-21 ± 4) . (7.3) 
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Let us finally discuss the contributions to ajf (Z, q) which are the ones relevant to the topic of 
this paper. Here, it is convenient to treat the contributions from the three generations separately. 
The contribution from the third generation can be calculated in a straightforward way with the 
result [El |2l 



EW(2) 



(r, t, b) 



G F m* a 
s/2 8tt 2 tt 
G F ra 
71 8^ 



-3 log 



2 



log 



Ml 
2 

6 



mt 



8 mj 
3 l0g M| 



•0 



M| m\ 



^3-7" x (-30.6). 



(7.4) 



In fact the terms of O log and ("^fj have also been calculated in ref. Q. There are 
in principle QCD perturbative corrections to this estimate, which have not been calculated, but the 
result in Eq. (7.4) is good enough for the accuracy required at present. 

As emphasized in ref. p2| , an appropriate QCD calculation when the quark in the loop of Fig. 1 is 
a light quark should take into account the dominant effects of spontaneous chiral symmetry breaking. 
Since this involves the u and d quarks, as well as the second generation s quark, we lump together 
the contributions from the first and second generation with the result, 



<^ 2 >(e, M , U ,d, 



s,c 



G F m- 
71 8^ 



-3 log z 



m 
M 
m 



-Slog^ 



M 



11 



4 log — f - — + -ir- 



7) If. 



6 9 



4, Ml 



■O 



Ml 



log 



Ml 



rr?„ a 



%^-x(-28.5±l. 
a/2 8tt 2 - 



+4.57 ± 1.80 + 0.04 ±0.02 I (7.5) 

(7.6) 



where the first line in Eq. (7.5) shows the result from the e-loop and the second line the result from the 
/i-loop and the c-quark, which is treated as a heavy quark (m c = 1.5 GeV). The first term in brackets 
in the third line is the one induced by the anomalous term in the hadronic < VAV > Green's function. 
The second term is the one induced by the leading effects of explicit chiral symmetry breaking in the 
non-anomalous longitudinal component of the < VAV > Green's function discussed in section 6 (the 
result in Eq. (f3J^)). The third term is the one induced by the transverse component of the < VAV > 
Green's function, evaluated in the pres ence of the light quark masses and in the HA with 3 poles to 
large-iV c QCD (the result in Eq. fl6.13| )). 



We want to stress that our result in Eq. (7.5) for the contribution from the first and second 
generations of quarks and leptons is conceptually very different to the corresponding one proposed in 
ref. [El, 



af W{2) {e,fJ<,u,d,s,c) = 



ml a 



G F 

V2S^ 



-3 log 



Ml 



01 M 7 

-31og^ 



4 log 



Ml 



, M z 



-i- 



Ml . Ml 11 
41og— |-log— f - — 



G F m 
71 8^ 



., x(-31.9). 



— 7T - 6 



(7.7) 



where the light quarks are, arbitrarily, treated the same way as heavy quarks, with m u = m e i = 
0.3 GeV, and m s = 0.5 GeV. Although, numerically, the two results turn out to be not too different, 
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the result in Eq. ( |7.7| ) follows from a hadronic model which is in contradiction with basic properties 
of QCD. The constituent quark model used to derive this result violates the current algebra Ward 
identity (A.l) derived in Appendix A. Furthermore, as discussed in the text, the model does not 
reproduce the QCD short-distance behaviour for the underlying < VAV > Green's function. These 
facts are at the origin of the spurious cancellation of the log Mz terms in Eq. ( [7.7| ) . A more detailed 
discussion of this issue can be found in Appendix D. 

Putting together the numerical results in Eqs. ( |7.2| ), ( |7.3D , (7.4) with the new result in Eq. ( |7.6| ), 
we finally obtain the value 



EW _ TO M 



y/2 8tt 2 



^ + i(l-4sin 2 M 2 



- (159 ± 4) 

7T 



(15.2 ± 0.1) x 10 



-10 



which shows that the two-loop correction represents indeed a reduction of the one-loop result by an 
amount of 22% . The final error here does not include higher order electroweak estimates Q. 
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APPENDIX 

A Ward Identities and Gauge Invar iance 



The hadronic Green's function in Eq. (2.1) is a combination of QCD three-point functions involving 
the vector and axial vector currents 

K W = (Vy 7^) (*) . and (x) = (V y 7^) (x) , 

of the algebra of currents, where A a denotes the flavour Gell-Mann matrices defined after Eq. ( |4.6| ). 
More precisely, with 

^V(?!,g 2 ) - i J d 4 x 1 d i x 2 e l ^ + ^ < n\T{v;-( Xl )vr(x2)A a a (o)}\n >, 

where 



the relation is the following 



A.l Current Algebra Ward Identities 

For a = 0, 3, 8, one has the following Ward identities |^| 

(qi+q 2 ) a T^ [a (q u q 2 ) = - J^L C a e^ p<x q^ 

+ J d 4 Xl d 4 x 2 e^-^+i*-^ < n | t{v™( Xi )v™(x 2 )d%(o)} I n > 

+ i tr (y) X / dix ^ x ^ l(qi ' X1+q2 ' X2) <"l T {^r(^i)^° m (^)(GG)(0)}|r!>,(A.l) 
with C 3 = 1, C 8 = I/a/3, C° = 2^2/3. Furthermore, 



D° = d abc M b P c , with P a = Vi75 — ^ : 



A-' 
2 

and 



1 [2 

M 3 =m u -m d , M S = -j={rn u + m d - 2m s ) , M Q = \J -(m u + m d + m s ) . 

Also 

8 

(GG) = a s e pvpa £ G (A)^ G (^W ) 
with G^ A > lil ' the gluon field strength tensor. Owing to the fact that 



K;?2> y ^irf 4 x 2 e^ 1 ^ +92 ^) <fi|T{y;™(x 1 )K om (^ 2 )(GG)(0)}|fi> = 0, 

f d i x 1 d i x 2 e l ^+^ <n\T{V-(x 1 )Vr(x 2 )D a 5 (0)}\n> = 0, 



10 Recall that we are using the following conventions: eoi23 = +1 and 75 = i7°7 1 7 2 7 3 , so that ^(7570717273) = 4i. 
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the above Ward identity can be solved as follows: 



(<7l + Q2)a f N c 



(<Zl+<7 2 ) 2 I 12^ 2 



C" 



+ ^2 tr (y) / d^d^e^-^-^ <n\T{v-( Xl )vr(x2)(GG)(0)}\n> 

+ J d 4 Xl d 4 x 2 e %^+^) < fi I T{V^{ Xl )Vr(x 2 )D? o {0)} \ Q > } 

+ ^; Q (<71,92), (A.2) 

where the tensor T a ^- a (qi, q 2 ) is now completely transverse, 

tii;q%\ {qi + q2) a }T^. a ( qu q 2 ) = o. 



The first term in Eq. ( A.2) generates the anomaly term discussed in the text. The second term is the 
one induced by the axial flavour-singlet U(l) anomaly. This term generates subleading contributions 
in the l/iV c -expansion with respect to the other terms; however, because of its flavour singlet nature, 
the numerical size of its contribution may not necessarily obey the expected l/iV c -power suppression. 
However, in our case, the source of the hadronic axial current is the one induced by the Z-boson 
coupling in the Standard Electroweak Theory, which has a vanishing flavour trace for each generation. 
Therefore, once the total sum of light and heavy quarks is taken into account, this term does not 
contribute to the muon anomaly. 



The third term in Eq. (A.2) is the one induced by the explicit breaking of chiral symmetry by the 
light quark masses in the QCD Lagrangian. We keep it, because it will be pertinent in the forthcoming 
discussions of electroweak gauge invariance. Notice, for further reference, that 

L>5 H ^= — Dg = m u uij 5 u - m d dij 5 d - m s si'y 5 s . 

V3 v6 



A.2 Electroweak Gauge Invariance 

Due to the conservation of the electromagnetic current, the gauge dependent terms in the photon 
propagator in the diagrams of Fig. 1 trigger naive Ward identities on the hadronic 7 — 7 — Z vertex, 
and therefore we can take the Feynman gauge for the photon propagator without loss of generality, 
as we have done in the calculation reported in the text. By contrast, the gauge dependence of the 
Z-propagator requires some discussion. In the Electroweak Theory, it is convenient to rewrite the 
Ward identity discussed in the previous subsection in terms of the coupling of the unphysical Higgs 
field <&°(:r). More precisely, in the notation used in the text, 

(<? - kfW^piq, k) = -i^C^ paf} q a k? - 2iM z W^(q, k) . (A.3) 

where 

C - Va e 2 with a g = +l/2, c q = +2/3 , for q = u,c,t\ 

C had -^a 9 e ? , with a 9 = _i/ 2 , e q = -l/3, for q = d,s,bj> {AA) 

and 

Wj$ 0) (q,k) = J d 4 xe^J tfive- a *{n\T{VZ"{x)V-(3,)J*°{Q)}\Sl), 

with 

T*° 1 \ ^ m q -■ 

J = 2 2^77-™. 
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This means that in full generality, in the Electroweak Theory, 

W^q, fc) = -iA (q 2 , fc 2 , (q - fc) 2 ) 7^|fc e Mpa/J q a k? + W£?{q, k) (A.5) 

with 

A (q 2 , k 2 , (q - fc) 2 ) = ^C had - 2M Z H* (q 2 ,k 2 , (q - fc) 2 ) , (A.6) 

where we have set 

W^ 0) (q, fc) = e, pa , (q 2 , fc 2 , (q - fc) 2 ) . 

Now, in the calculation corresponding to Fig. 1, the most general form of the Z-propagator, in a 
renormalizable linear gauge of the 't Hooft class, is 



(q-ky(q-ky' 
{q-kf-i z Ml 



-r'(q-k) 2 -{q-kY(q-kY'} - &M 2 



(q-kY( q -kY 
Ml 



[{q-kf-M 2 ] [(q-kY-tzM 2 ] 



(A.7) 



(A. 



The unitary gauge, which was used e.g. in ref. ]22[ ], is recovered in the limit £z ~~ > °°7 which has to be 
taken from the start. In the text we have done the calculation in the Feynman gauge £z = 1- Clearly, 
there is no contribution to the muon anomaly from the contraction of the gauge dependent piece in 
the Z-propagator in (A.7) with the fully transverse W flvp (q,k) hadronic tensor. The contribution 
from the non-transverse term in Eq. ( |A.5| ) however, is gauge dependent. It selects the second term in 
the iv-propagator in (A. 8), via the contraction 



<2 



(q - fc) 2 (q - fc) 2 - izM\ 



(g ~ W 



As shown in Eq. (A6) the non-transverse term in Eq. (|A.5| ) has two pieces. Let us first discuss the 
one coming from the anomaly term. For a specific fermion /, the piece coming from the anomaly 
results in the expression (in the £z gauge) 



F 2 (0) 



(/) 



G F m a 



In 



izM 2 



while in the unitarity gauge, the corresponding result is divergent and requires regularization. In 



dimensional regularization, with A^s 



■ ■ 



l 

ci-4 



;iog47r + 2 + r'(l)) 



one finds 



F 2 {0) 



(/) 





i" 




+ 2 



Recall that in these expressions, when the fermion / is a quark q, a 9 and e q have been defined in 
Eq. (A. 4). When / is a lepton I, &i — —1/2 and e; = —1. In particular, for f = u,d,s , 



u,d,s 



4JV„ 



1 



(-1) 



Nc 4 
3 3' 



and in the Feynman gauge where £,z = 1, we recover the result of Eq. (3.5). Notice also that, in the 
Standard Theory, the £z~gauge dependence, or the log {^r^j — 327t 2 Ams dependence in the unitary 
gauge, cancel when quarks and leptons are taken together generation by generation. 
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Next, we discuss the gauge dependent piece coming from the term proportional to 7i* (q 2 , k 2 , (q — ky 
in Eq. (A.6), i.e. the term 



V^ z) (q-k)xH 



$ o (q - k) v > $ o -i£ z (q - k) v 

— n 



(q - ky 



{q-ky-tzM 2 (q-k) 2 



This term, when acting on the leptonic line (see Eq. (1.3)) triggers a trivial Ward identity: 



(q - k) v 

-2m jLl u(p / ) 



7 M -» 7 — — 7^75 + 7^75 —, 1 — 7 M 



- 9 - m„ 



4+ d-m. 



*(p) = 



u 1 1 a 

*)'- d - m,, P+ - m,. 



u(p) 



(A.9) 



The resulting expression: 
-Hz 



(q - k) 2 - i z Ml (q - k) 



x 2m fl u(p') 



„ % 1 
T u_ 4 — 1—75 + 75 



r 



u(p) 



when added together with the gauge dependent term from the unphysical Higgs contribution to the 
muon anomaly, 



M 2 [{q-k) 2 -i z M 2 



x 2m ll u(p') 



„ 1 1 i j. 

V— 9— m u P+ d - m,. 



u[p) 



results in an overall gauge independent contribution, proportional to 



M 2 {q-ky 



x 2m ll u(p') 



u(p). (A.10) 



In the unitary gauge, there are no unphysical Higgs couplings and the piece coming from the term 
proportional to in Eq. (A.6) generates the factor 



H 



fJ 



viz' _ {q-k) v (q-k)" 
Ml 



(q - k) 2 - M 2 Z 



(g - k) v , _ $0 i 

(q - k) 2 n M 2 (q - ky- 



:(q-ky, 



which, after contraction with the leptonic line, produces an expression identical to the one in Eq. (A. 10) 



B Technical Remarks on the OPE in the Chiral Limit 

As already mentioned in the text, the relevant OPE is the one of the product of the two currents 



V^ m (a;)A^ c (2/) in Eq. (2.1) when x — > y, because these are the currents where the hard virtual momen- 



tum q can get through. There are two types of contributions to discuss in the chiral limit. 

B.l Contribution from Soft Gluons 

These are the contributions generated by the term 

lim f d A ze^ z T{V-(z)A-(0)}^ 

q^oo J I" 

: ^(0)7^ Q iS iS {l)lu^Qf q J (0) ■ + ■ q 3 (0) lul5 Q^ iS ji {q)^Q <f (0) : (B.l) 
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in the perturbative Wick-expansion denote quark color indices) when the short-distance expan- 
sion of the quark propagator ^ is inserted: 



(B.2) 



iS t3 (q) = d k ze t ' !Z iS l: >(z,0) = 5 v ~ 



i 



~gX [flT&lfff- 4 D a G% - q^D, q«G%^ - Mq^D^G%^^\ (B.3) 
o q 



where ( A is a color index in the adjoint representation) 

Gq/3 = 9 X! G^ (0) ) G^ = -(-ajipaG 



pa 

ij 



The terms of lowest dimension in the short-distance expansion of the quark propagator in Eq. (B.2) 
which can contribute (a priori) to a tensor structure like the one in Eq. (4.2) are those with two 
uncontracted powers of g-momenta in the numerator i.e., the last three terms in Eq. ( |B.3j ). However, 
these terms, as far as their Dirac structure is concerned, are all proportional to one gamma matrix or 
to 7^75, (the last term.) Therefore, when inserted in the r.h.s. of Eq. (B.l), and taking into account 
that 



plus the identity 



QQf =Qf Q = idiag(2,l,l), 



IpIpIu = 9pplv + 9 P ulp - 9pul P - itppuS 7 75 



(B.4) 



one finds that these terms cannot produce an antisymmetric tensor in the /x and v indices. We then 
conclude that there is no contribution from soft-gluons with the tensor structure of Eq. (4.2) with a 
power p < 3. 



B.2 Contribution from Hard Gluons 



Expanding the T-product in Eq. (4T) in perturbation theory to 0{gf) generates the following operator 
U» v {q) = ^9 2 s Um / d 4 z e^ z J d 4 x J d 4 y x 

r j: Ux)l a ^fGi A Xx)^ 3 {x) : : Mvh^G^ MMv) ■ 
: q(z)^Qq(z) : : qiOh^Q^kiO) } ■ 



The contraction of the two gluon fields generates four possible classes of hard topologies illustrated in 
Fig. 3 above. Using the identity 

and the fact that the second term in this identity generates subleading terms in the l/iV c -expansion, 
the contributions from the four hard topologies can be written as follows: 

"See e.g. Eq. (2.34) in ref. Notice that iS ji (q) = J d i ze- iq ' z iS : > i (0, z) . 
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q— s-oo 



q— >oo 



•« 2 2 1 
z! g u 
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« 2 2 1 
— i — g„ 
2r s g 6 


g af3 


{ft (S^a -+ 


-A lia )Qq i 
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+ A l , l3 )T 3l5 q t }{0) 


i 2 2 1 
2r s q 6 


g aP 


{q l (S^ - 


A m )Qft 


ft 


{Su0 


+ A, (3 )T 3 7 5 <f}(0) 


i 2 2 1 
2r s g 6 


g af3 


{ft {S m 4 


~A^ a )Qq3 


q j 


(Suf3 


-Aup)T 3l5 ^}{0) 



u${q) 

K 2 J(q) 
U$(q) 
U${q) 

where S^ a and A pa denote the tensors: 

S>a = q^la + q a lu, ~ 9u,a 4, and A pa = iq" 1 e maa ^ a 75 . 
Only the antisymmetric tensors contribute in the total sum, with the result 

zV(?)Uoc = ~^9 2 s ^g a0 ^) 2 q v q x e maa e vX P (q 4 l°l 5 Qft)(q j 7 P T 3 ft)(O) . 
Fierzing the resulting four-quark operator, using the relevant terms of the identity 

(SyWH^yV) = -\e afSpa [m{^a pa q) + {qa pa ft{ft)\ + ■ • • , 
and contracting two of the epsilon tensors (using the correct sign!) one finally gets the result 

U»v{q) = i [qs(-^u a pq a - q^p, va &q a \ ( -2tt 2 — h — ; 

with O p5 (0) the tensor 



0^(0) 



\ (ua^u) {uu) + \ (da^d) (dd) + \ (sa^s) {ss) 



(0). 



C The OPE beyond the Chiral Limit 

In the presence of quark masses, there are two types of terms that one has to consider in the evaluation 
of the function in Eq. (4.1). (Recall that the relevant terms in the OPE are those with the explicit 
tensor structure shown in Eq. ( |4.2j ).) There are terms of 0{m) from the expansion of the quark 
propagator: 

1 111^,,,, 
+ ~j m ~j + C( TO ) , 



4-m 4 4 4 

and terms from the Taylor expansion of the quark field 

q{x) = q{0) + X a D a q{0) + ■ ■ ■ . 
The first type expansion leads to a potential contribution: 

ZV(<?)| m = 1 ^2 m° a0 MQQfq{O) , (C.l) 

where M. denotes the quark mass matrix 

M = dia,g{m u ,m d ,m s ) . (C.2) 
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This contribution, however, does not have the two powers of g-momentum required by the tensor 
structure in Eq. (4.2). 

The Taylor expansion from the quark field leads to the operator 

U, v (q)\ q = D p q(0) 1 , Yq2 ~l iqP lvl5 QQfq{Q) - g(0) 7 ,75 ^ ~J ^ IhQTQ D p q(0) ■ (C.3) 



Using the Dirac equation 

(iD^ p -M)q = 0, 
the terms proportional to j p in Eq. (|C.3|) lead to the operator 



D^mi^QQfqiQ) - g(0)7,75<9Qi 3) ^g(0) 



(C.4) 



(C.5) 



where we have used the fact that terms ~ a a P cancel with those from the quark-propagator contri- 
bution in Eq. (C.l). The other terms in Eq. ( |C.3[ ) are proportional to ^ , and lead to the following 
combination of operators 

= - J* { iq^uasq a (D f3 q(o)i s QQL } m + m^QQ^D^qio)) 

+q f3 (q^ s - g plJ q 5 + qj^) (d p q{0) lslb QQf q{0) ~ q{Q)lS^QQf D q{Q)) } .(C.6) 



When inserting the operators U pv (q)\^ p and U pv {q)\^ in the T-product in the r.h.s. of Eq. (2.1) 



there appear the following two-point functions 

sP' p (k) = J d^e-* k -y(n\T { (zAV QQ [ l\ + qi 5 QQ { L 3) D^) (o)^T(y)} |n> . (c.7) 

and 

A^p(fc) = J d*ye-* k -y(n\T { (d p qjsToQQ^Q ~ qislsQQ^ D p q) (0)V™(y)} |fi) . (C.8) 

Lorentz covariance and gauge invariance of the electromagnetic current restrict these two-point func- 
tions to the following decomposition p] 



and 



Zpsp(k) = -(gspk 2 - k s k p ) fc /3 riy/y(fc 2 ) 



A /3 5p(fc) = iei3p8 a k a Tl A , v (k 2 ) . 



(C.9) 
(CIO) 



We can now obtain the asymptotic behaviour of the Green's function W llvp {q, k) defined in Eq. (2.1) 
in terms of the invariant functions Hyv(k 2 ) and IlA'v(k 2 ), with the result 



lim W llvp {q,k) 



2i 



y^£ l ipvak a Tl-Aiv{k 2 



2i 

~7t4 {(9M e /3pfQ + Qvtppua) q k a Ii A , v (k 2 ) 
+e liav fi {p <3 pk 2 - k S k fi ) q a (q ■ k)n v , v (k 2 )} 



(C.ll) 



12 There is in fact a second possible term in E.pg p (k) of the type — (o o„fc 2 — kpk p ) kg Tl v iy(k 2 ), which however vanishes 
as can be seen by contracting the /3 and p indices in the r.h.s of Eq. (0.7), using translation invariance and conservation 
ofV„ em . 
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The next step consists in extracting from the expression in Eq. (C.ll) the transverse component: 
lim^oo W^p a (q, k), defined in Eq. (2.6). Taking the derivative with respect to the fc-momentum, as 
defined in Eq. (O), we finally get 



2i 



lim W livpa {q) = --tt7 [q P e^aaq a - q<j£^ ap q a ] hU'v(O) . 

It follows then that the asymptotic behaviour of the function W(Q 2 ), in the presence of quark masses, 
is given by 

lim W{Q 2 ) =--^IW(0). 



What do we know a bout 11^4/^(0)? It turns out that the invariant functions H A'v jk 2 ), Hyv(k 2 ) 
defined in Eqs. (|C.9| ), ( |C.10 ) and the invariant function Tlvrik 2 ) defined in Eq. ( |4.6| ) are not inde- 
pendent. This can be seen starting with the identity 

iqcTapMQQ^ q = - ^D p qj p a a pQQ^\ - qa a p'-fpQQ ( L ) Di'q^) , 



which follows from the Dirac equation ( p.4[ ) , and using the Dirac algebra identity in Eq. ( BT ) which 
brings in the relevant operators. The resulting relation, to lowest order in the quark masses, is 

2tr (MQQfQ) n VT (fc 2 ) - \ k 2 U v , v {k 2 ) + IW(fc 2 ) , 

which at k 2 = reduces to 



rU'v(O) = g(4m„ ~m d ~ ra^II^O) 



D On the presence of log M z terms in a^ w( - 2 \e, fx, u, d, s, c) 

We have noticed in section 7.2 that for the two first generations, the log Mz terms do not cancel 
in ap^ W ^ (e, u, d, s, c), see Eq. (7.5). This result might appear as puzzling, since in the case of 



free quarks, which encompasses the constituent quark model, the cancellation of the log Mz terms in 
a^ W ^ (e, fi, u, d, s, c) occurs, as observed in ref. [^3), and as shown in Eq. ( [7.7[ ). However, this latter 
situation is conceptually very different from the one described by confined light (or even massless) 
quarks. This difference, that needs perhaps to be explained more in detail, is at the origin of the 
non-cancellation of the log Mz contributions in Eq. (7.£). 



In order to proceed, let us start from Eq. (5.7), that we write in a slightly different, but more 
convenient, form, 



^xqmO) 



N e 8 M 2 



12tt 2 3 
dt 



^ 9A w x 1 dy 



1 V dx Xy ~ - ^ 
X l + j^y{l-y) 



Q 



t + Q 2 — ie 7T 



lmw x q M (t). 



(D.l) 



The second, i.e. a dispersive form, is obtained upon performing the integration over the Feynman 
parameter x in the expression on the first line, followed by an obvious change of variable in the 
remaining integration over y, with the outcome 



-lmw x q M (t) 



XMl 
2 t 2 



1 



(D.2) 
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For the case of a massless free fermion, one obtains 



— N 2 M 2 

Jim oWxQM (z) = Y ^- 9AW ^- , (D.3) 



since, as is well known, 



1 M 2 1 1 

lim - - , 6{t - AMI) = — 7 5{t) . (D.4) 

m q ^o 2 t 2 r Tm% v q > 4 v ; v 7 



t 

Thus, in the massless limit, the function u^qmM reduces to a simple pole. This singularity at Q 2 = 
also provides the 1/Q 2 behaviour at large Q 2 , which in turn will lead to the cancellation of the log Mz 
terms. Notice also that this asymptotic 1/Q 2 behaviour, and the ensuing cancellation of log Mz terms, 
is preserved for a massive fermion, as shown by Eq. ( f5.7| ) or by the dispersive representation above. 

The real question we now have to address is how much of all this can be taken over to QCD. Let 
us first discuss the chiral limit. One immediate and important observation one has to make in the 
case of QCD is that the pole we have just discussed in the case of a massless free fermion no longer 
occurs. Indeed, one may evaluate the relevant three-point function W l _u, p {qi,q2) within the effective 
low-energy chiral lagrangian. The lowest order contribution starts at order 0(p 4 ) and is given by 
the anomaly pole, as discussed in ref. p2]] . At NLO, 0(p 6 ) in this case, we have a counterterm 
contribution ^|, which leads to a non vanishing, but constant, value of w(Q 2 ) Thus, confinement 
drastically modifies the long distance behaviour of W(Q 2 ) as compared to the free fermion case, not 
really a surprise. It now remains to discuss the short distance aspects. We may write the function 
w(Q 2 ) in the dispersive form given before for the free fermion case, 

f°° II 

w(Q 2 ) = dt—— 2 --lmw(t), (D.5) 

J-ao t + Q z ~ie 7T 

but where now the perturbative QCD part describes the absorbtive part above a certain threshold 
so, which has to be sufficiently large so that perturbative QCD can be applied. Within the large- Nc 
framework, we may thus write 

-lmw{t) = -lmw Res {t) + -lmw pQCD {t). (D.6) 

TT IT 7T 

The non perturbative piece from the narrow width resonant states writes 

-lmw Res (t) = V M 2 a R S(t - M 2 ) , (D.7) 

7T * — ' 

R 

and Mr the mass of the R narrow state, M R < sq. The perturbative contribution for a massive light 
quark with a current algebra mass m q , to lowest order in perturbative QCD (free quarks, <?a = 1), is 

I Im ^co (i) _ _ 1 £ Ht _ „) (^M°) . (D.8, 



This gives then 

w (z) = w Res (z) + w pQCD (z) , (D.9) 

13 We stay of within the large- Nc framework, hence meson loops can be ignored. 

14 In terms of the C(p 6 ) chiral Lagrangian constructed for the odd intrinsic parity sector in the relevant low-energy 
constant is C^. 



27 



with 



w Res (z) 
w pQCD (z) 



E 



Z + Zr 



zr 



Ml 



1^ 
4 Q- 



Am 2 , 
so 



log 1 + — 

Q l \ so 



0{m\) 



Nr 



12tt 2 3* 



-g A \(l + 0(a s )) 



We see here the crucial difference with the constituent quark model, where Sq would coincide with 
4Mq, thus providing the same leading 1/Q 2 short distance behaviour as for a massless free quark. In 
QCD, the onset of the continuum occurs at a scale so that is much larger than the scale set by the 
light quark masses, and for sq 3> 4m 2 , one finds 

lim zw pQCD (z) = G(ml). 

z — >oo 

Thus, in the chiral limit, or if one considers only first order explicit chiral symmetry breaking effects, 
the short distance behaviour of w(z) in QCD reads as given in Eq. ( |6.10 ), 



r . / \ 4 (4m„ - m d - m s ){ipif)) 
lim wiz) = -. — — 
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2 q s 

— 7T + ■ 

9 7T ne 



The fact that the leading short distance fall off of w(z) in QCD is 1/z 3 in the chiral limit, or 1/z 2 
for massive light quark if only first order order explicit chiral symmetry breaking effects are retained, 
explains why the cancellation of log Mz terms no longer occurs in QCD. 

We conclude therefore that the cancellation of the logAfz terms does not occur in QCD in the 
chiral limit, and by stressing once again very strongly that constituent quark models do not, in general, 
provide an adequate description of fundamental properties of QCD. 
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